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MODELING  THE  PRESSURE-DILATATION  CORRELATION  1 


S.  Sarkar 

Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA  23665 

ABSTRACT 

It  is  generally  accepted  that  the  pressure-dilatation,  which  is  an  additional  compressibility 
term  in  the  turbulence  transport  equations,  may  be  important  for  high-speed  flows.  Recent 
direct  simulations  of  homogeneous  shear  turbulence  have  given  concrete  evidence  that  the 
pressure-dilatation  is  important  insofar  that  it  contributes  to  the  reduced  growth  of  turbulent 
kinetic  energy  due  to  compressibility  effects.  The  present  work  addresses  the  problem  of 
modeling  the  pressure-dilatation.  We  first  isolate  a  component  of  the  pressure-dilatation 
which  exhibits  temporal  oscillations  and.  using  direct  numerical  simulations  of  homogeneous 
shear  turbulence  and  isotropic  turbulence,  show  that  it  has  a  negligible  contribution  to 
the  evolution  of  turbulent  kinetic  energy.  Then,  an  analysis  for  the  case  of  homogeneous 
turbulence  is  performed  to  obtain  a  model  for  the  non-oscillatory  pressure-dilatation.  This 
model  algebraically  relates  the  pressure-dilatation  to  quantities  traditionally  obtained  in 
incompressible  turbulence  closures.  The  model  is  validated  by  direct  comparison  with  the 
pressure-dilatation  data  obtained  from  the  simulations. 


1  This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


1  Introduction 


The  pressure-dilatation  appears  as  an  explicit  compressibility  term  in  the  equations  for  mean 
temperature  and  turbulent  kinetic  energy.  It  has  been  generally  recognized  that  this  term 
may  be  large  in  high-speed  flows  and  therefore  requires  consideration  in  compressible  turbu¬ 
lence  closures.  In  Sarkar,  Erlebacher,  Hussaini  and  Kreiss1,  it  was  found  that  the  pressure 
dilatation  was  smaller  than  the  compressible  dissipation  in  direct  simulations  of  isotropic 
compressible  turbulence  which  suggested  that  the  presure-dilatation  can  be  absorbed  into 
the  model  for  the  compressible  dissipation  derived  therein.  However,  subsequent  direct  sim¬ 
ulations  of  homogeneous  shear  flow  by  Blaisdell  et  al.2  and  Sarkar,  Erlebacher  and  Hussaini3 
showed  that  the  pressure-dilatation  is  comparable  to  the  compressible  dissipation  and  con¬ 
tributes  to  the  reduced  growth  of  turbulent  kinetic  energy  induced  by  compressibility.  This 
has  motivated  a  revisit  to  the  issue  of  modeling  the  pressure-dilatation  in  the  present  paper. 
Recently,  Taulbee  and  VanOsdol4  and  Zeman5  have  also  considered  this  problem.  T^.ulbee 
and  VanOsdol  related  the  sum  of  the  pressure-dilatation  and  compressible  dissipation  to  a 
model  involving  the  density  variance  and  the  divergence  of  the  mean  velocity.  A  separate 
modeled  transport  equation  was  proposed  for  the  density  variance.  Zeman5  equated  the 
pressure-dilatation  to  the  time  rate  of  change  of  pressure-variance  and  then,  extending  the 
work  by  Sarkar  et  al.3  on  equilibration  of  the  compressible  pressure  variance  on  the  fast 
acoustic  time  scale,  proposed  a  transport  equation  for  the  pressure  variance.  The  objective 
of  the  present  work  is  to  deduce  a  model  for  the  pressure-dilatation  which  is  algebraically 
related  to  quantities  obtained  bv  incompressible  turbulence  closures,  in  contrast  to  the  pre¬ 
vious  models,  which  require  new  transport  equations. 

We  consider  the  behavior  of  the  pressure-dilatation  in  homogeneous  flows.  Direct  simu¬ 
lation  results  are  presented  and  compared  for  two  flows  -  homogeneous  shear  turbulence  and 
decaying  isotropic  turbulence.  The  DNS  results  show  that  the  contribution  of  the  pressure- 
dilatation  to  the  turbulent  kinetic  energy  evolution  is  more  important  in  the  case  of  shear 
turbulence  than  in  decaying  isotropic  turbulence.  We  also  find  the  somewhat  surprising  re¬ 
sult  that  the  major  contributor  to  the  pressure-dilatation  is  not  the  compressible  pressure, 
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but  is  the  incompressible  pressure  associated  with  the  solenoidal  velocity.  The  governing 
equations  are  analyzed  for  the  case  of  homogeneous  flow  to  obtain  a  formal  expression  for 
the  pressure-dilatation  from  which  a  model  is  deduced  using  scaling  arguments. 

2  DNS  results 

A  spectral  collocation  method  along  with  a  third-order  accurate  Runge-Kutta  time  advance¬ 
ment  was  used  to  solve  the  compressible  Navier-Stokes  equations.  Details  of  the  calculations 
are  available  in  Sarkar,  Erlebacher  and  Hussaini3  for  homogeneous  shear  flow,  and  in  Er- 
lebacher,  Hussaini,  Kreiss  and  Sarkar6  for  decaying  isotropic  turbulence.  A  uniform  963 
grid  was  used  for  discretizing  the  flow  domain.  Table  1  shows  some  of  the  parameters  of  the 
three  shear  flow  cases  (Sl,S2  and  S3),  and  three  decaying  isotropic  cases  (D1,D2,  and  D3).  In 
Tables  1  and  2,  5  and  r  denote  the  shear  rate  and  kinematic  viscosity  used  in  th^  (nondimen- 
sional)  compressible  Navier-Stokes  equations,  while  R\  0  and  Mt  0  denote  the  initial  values 
of  the  Taylor  microscale  Reynolds  number  and  turbulent  Mach  number  respectively.  Note 
that  Re\  =  qX/u  where  q  =  \Ju'xu[  and  A  =  q/Ju>[oj[  ,  while  Mt  =  q/c  where  c  =  yyRT  is 
the  mean  speed  of  sound.  In  order  to  minimize  the  introduction  of  compressible  effects  due 
to  initial  conditions,  all  the  cases  start  with  incompressible  data,  that  is,  the  velocity  field 
is  divergence- free  ( d "  =  V.u'  =  0),  the  pressure  field  satisfies  the  usual  Poisson  equation  for 
incompressible  flows,  and  density  fluctuations  p'  =  0.  The  initial  temperature  is  calculated 
from  the  equation  of  state  using  the  known  pressure  and  density  fields. 

The  homogeneous  flows  considered  here  have  temporally  evolving  turbulence  statistics. 
Fig.  1  shows  the  evolution  of  the  pressure-dilatation  p'd'  as  a  function  of  normalized  time 
St  for  the  two  homogeneous  shear  cases  SI  and  S2.  Two  trends  in  the  behavior  of  p'd' 
are  evident  from  Figs,  la-b;  first,  the  pressure-dilatation  develops  pronounced  oscillations 
in  time  and  second,  it  is  more  negative  than  positive.  These  trends  are  common  to  all  the 
simulations  of  homogeneous  shear  flow  that  we  have  performed.  Such  oscillations  in  statisti¬ 
cal  quantities  arc  neither  expected  nor  encountered  in  incompressible  flows  and  complicate 
issues  in  compressible  turbulence  modeling.  By  numerical  experiments,  it  was  found  that 
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the  nominal  time  period  of  the  oscillations  decreased  approximately  linearly  with  the  speed 


of  sound.  This  suggested  that  one  could  isolate  the  oscillatory  part  of  p'd'  by  decomposing 
the  fluctuating  pressure  p'  into  the  sum  of  an  incompressible  part  p1'  and  a  compressible 
part  pc' .  The  component  p1'  is  associated  with  the  incompressible  velocity  field  u1  which  is 
divergence-free  (V.u7  =  0)  and  satisfies  the  usual  Poisson  equation 


VV'  =  -2 P^y'n^  - 

and  the  remainder  p'  —  p1'  is  the  compressible  pressure  pc> .  Since 

p  =  p1'  +  pC' 


(1) 


we  have 


p'd,'  =  p1' d'  pc‘  d' 


Fig.  2a  shows  the  evolution  of  pc'd'  and  p1' d'  for  case  Si.  The  oscillations  are  substantial 


only  for  pc' d' ,  and  furthermore,  the  peaks  and  valleys  in  the  evolution  of  pc' d'  in  Fig.  2a 
seem  to  be  much  more  symmetric  around  the  origin  than  those  in  p'd'  in  Fig.  la. 


One  of  the  important  effects  of  p'd'  is  its  influence  on  the  budget  of  the  turbulent  kinetic 
energy  K .  The  equation  for  K  in  homogeneous  turbulence  is 

d 


=  ~pP  ~  pc,-  pcc  +  p'd' 


(2) 


where  V  =  — is  the  production,  e,  the  solenoidal  dissipation,  tc  the  compressible 
dissipation,  and  ptd'  the  pressure-dilatation.  The  overbar  denotes  a  conventional  Reynolds 
average,  while  the  overtilde  denotes  a  Favre  average.  A  single  superscript '  represents  fluctua¬ 
tions  with  respect  to  the  Reynolds  average,  while  a  double  superscript  "  signifies  fluctuations 
with  respect  to  the  Favre  average.  In  order  to  gauge  the  relative  importance  of  the  two  com¬ 


ponents  pc'd'  and  pJ  d'  of  the  pressure-dilatation  in  the  evolution  of  the  turbulent  kinetic 


energy,  we  plot  the  integrals  /  p1' d'  dt  and  /  pc‘ d'  dt  in  Fig.  2b.  Fig.  2b  shows  that  the  inte¬ 


grated  contributor  of  p1' d'  is  about  an  wu<u  of  magnitude  Ictiger  than  ti.u*  of  p^'cf.  Thus, 


even  though  the  extrema  of  pc'd'  are  larger  than  p1' d\  the  time-integrated  contribution  of 
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these  oscillations  in  pc‘ d'  is  negligible  compared  to  the  contribution  of  p1' d’ .  Fig.  3  and 
Fig.  4,  which  contrast  the  behavior  of  p1' d'  and  pc'd'  for  two  other  homogeneous  shear  flow 
cases  S‘2  and  S3,  confirm  the  aforementioned  trends  observed  in  case  SI.  (See  Table  1  for 
the  initial  conditions  of  cases  S2  and  S3).  It  seems  that  for  homogeneous  shear  flow,  p1' d' 
dominates  pc‘  d'  in  the  kinetic  energy  equation.  It  is  interesting  that  even  though  both  the 
dilatation  and  the  compressible  pressure  are  associated  with  the  same  hyperbolic  system  of 
equations  governing  the  compressible  mode,  the  correlation  of  the  dilatation  with  the  com¬ 
pressible  pressure  plays  a  smaller  overall  role  than  the  correlation  of  the  dilatation  with  the 
incompressible  pressure. 

Figs.  5  and  6  show  DNS  results  on  the  pressure-dilatation  from  the  decaying  isotropic 
turbulence  cases  Dl  and  D2  (see  Table  2  for  initial  conditions).  Both  p1' d'  and  pc' d'  in 
Fig.  5a  and  6a  show  sharp  transients  initially  whereby  the  pressure  field  is  appropriately 
redistributed  into  incompressible  and  compressible  parts.  The  integrated  p1' d'  is  much  larger 
than  the  integrated  pc'd'  in  Figs.  5b  and  6b  just  as  in  the  homogeneous  shear  flow  case.  In 
Figs.  5a  and  6a,  the  term  pc'd'  becomes  approximately  zero  after  a  short  initial  transient. 
This  is  in  agreement  with  the  exact  solution  derived  in  Sarkar  et  al.1  for  the  linearized 
equations  for  the  compressible  mode  which  predicted  that  pc' d'  — ►  0  after  a  transient  on 
the  acoustic  time  scale.  The  DNS  reveal  an  important  difference  in  the  behavior  of  p/d! 
in  isotropic  turbulence  with  respect  to  homogeneous  shear.  The  term  p'd'  is  predominantly 
positive  in  the  case  of  decaying  isotropic  turbulence,  in  contrast  to  its  predominantly  negative 
behavior  in  homogeneous  shear  turbulence.  The  different  signs  of  the  pressure-dilatation  have 
been  explained  by  a  theoretical  consideration  of  the  equations  of  the  pressure  variance  and 
density  variance  by  Sarkar,  Erlebacher  and  Hussaini3. 

To  summarize,  the  DNS  results  show  that  the  pressure-dilatation  is  predominantly  neg¬ 
ative  and  has  pronounced  oscillations  in  homogeneous  shear  flow,  while  it  is  predominantly 
positive  in  decaying  isotropic  turbulence.  We  find  that  the  oscillations  in  p’d'  are  confined 
to  the  component  pc' d'  of  the  pressure-dilatation  associated  with  the  compressible  pressure. 
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and  furthermore  that,  because  of  self-cancellation,  the  contribution  of  this  oscillatory  com¬ 
ponent  to  the  development  of  the  kinetic  energy  is  negligible.  Therefore,  it  seems  that  only 
the  component  p1'  d'  of  the  pressure-dilatation  requires  modeling,  which  we  proceed  to  do  in 
the  next  section. 


3  Modeling  the  pressure-dilatation 


Consider  the  Poisson  equation  Eq.  (1)  for  the  incompressible  pressure.  As  is  often  done  for 
pressure-strain  modeling  (see  Lumley7,  Reynolds8)  for  incompressible  flows,  the  incompress¬ 
ible  pressure  can  be  split  into  a  rapid  part  pR'  which  reacts  instantaneously  to  a  change  in 
the  mean  velocity  gradient  and  a  slow  part  ps' .  Thus 


/  =  PR' 


+  P" 


(3) 


where 


\~/2R'  o—-  /  V 

V  V  =  -2 pumnu  n>m 


u) 


ana 


V72  S'  -  I'  /' 
V  p  =-pu  m  nu  n. 


(5) 


Eq.  (4)  can  be  exactly  solved  by  Fourier  transforms  for  homogeneous  turbulence  to  give 


R'  o  —  -  i  « i 

P  =  2lPjYum,nlln 


•_  /  a  / 


(6) 


Note  that  (f>  denotes  the  Fourier  transform  of  <f>.  From  Ecp  (6)  it  follows  that 

p**  =  2m!.»J^rE%jk  (-) 

Here  denotes  the  spectrum  of  the  mixed  Reynolds  stress  tensor  .  Eq.  (7)  is  an 

exact  expression  for  the  rapid  pressure-dilatation  pR' d'  and,  though  cumbersome,  can  be 
used  for  obtaining  a  simple  model.  By  inspection  of  Eq.  (7),  it  is  clear  that  because  of 
the  dependence  on  the  local  energy  spectrum  tensors,  a  transport  equation  is  required  for 


a  general  representation  of  p^'d'.  However,  we  will  attempt  to  obtain  an  algebraic  rather 


than  a  differential  model  for  compressibility  terms.  This  is  possible  if  there  are  equilibrium 
scalings  in  the  flow,  and  the  ensuing  model  will  be  useful  if  the  additional  compressibility 
correlations  do  not  dominate  the  incompressible  terms  in  a  given  transport  equation. 

Define  the  tensor  Amn  by 

=  j  (8) 

which  enables  the  following  compact  representation  of  Eq.  (7)  of  the  rapid  pressure-dilatation 

pH'd'  =  2puImnAmn  (9) 

The  simplest  dimensionally  consistent  form  for  Amn  which  has  the  correct  dependence 
on  u\'  and  uf'  is 

Amn  =  a2u£'u''  (10) 

where  a2  is  a  dimensionless  parameter  which  in  general  is  a  function  of  the  actual  shape 
of  the  energy  spectrum  tensors  for  the  incompressible  and  compressible  velocity.  We  have 
found  through  analysis  and  DNS  that  for  homogeneous  shear  turbulence 

uc=0{Mt)uI  (11) 

where  uc  and  u1  respectively  denote  the  L2  norms  of  the  compressible  and  incompressible 

velocity  fields.  In  Sarkar  et  aid,  we  had  shown  that  the  dilatational  velocity  has  a  fast  time 
scale  which  is  0(Mt)  times  the  solenoidal  velocity’s  time  scale,  and  therefore  the  correlation 
between  u'J  and  tt£  should  be  prorated  by  a  factor  of  Mt.  Using  Eq.  (11)  and  prorating 
the  mixed  correlation,  we  obtain  from  Eq.  (10), 

Amn-  ^MtU'mU'n  (12) 

Substituting  Eq.  (12)  into  Eq.  (9)  gives 

p^d'  =  ~a2M?V  (13) 

where  V  is  the  production  -w4X„«'n-  Eq.  (13)  may  be  used  as  a  model  for  the  rapid- 
pressure  dilatation  when  using  a  conventional  Reynolds-averaged  system  of  equations.  We 
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prefer  using  a  combination  of  Favre  and  Reynolds  averages  (eg.  see  Sarkar  and  Balakrishnan9 
for  .ompressible  flow  calculations,  and  for  such  a  sysiem  of  equations  we  propose  using 
Ecp  (13)  with  the  production  being  V  —  —  Since  p' /p  —  0(A/t2),  changing  from  the 

Reynolds-averaged  definition  to  the  Favre-averaged  definition  of  the  production  V  introduces 
terms  of  higher  order  in  Mt  which  may  be  neglected. 

Let  us  now  consider  the  remaining  part  of  the  pressure-dilatation,  the  slow  pressure- 
dilatation  ps'd'.  After  using  Fourier  transforms  to  solve  Eq.  (5)  for  the  slow  pressure,  we 
obtain  the  following  expression  for  the  slow  pressure-dilatation, 


P 


(1-0 


Here  <£*  denotes  the  complex  conjugate  of  the  Fourier  transform  <p.  An  order  of  magnitude 
analysis  of  the  r.h.s.  of  Eq.  (14)  gives 

/V  =  ?0(A)0(A)a/|  (15) 

u‘  l 

In  Eq.  (15),  /  denotes  the  integral  length  scale  of  the  turbulence,  and  the  last  factor  Mt  arises 
from  the  disparity  between  the  scales  of  the  incompressible  and  compressible  fields.  Using 
Eq.  (11)  for  the  scaling  of  uc ,  and  noting  that  t3  —  0(ul3 //),  we  obtain  the  model 

pS'd'  =  a3pc3M?  (16) 

Combining  Eqs.(13)  and  (16),  we  have  the  following  model  for  the  incompressible  pressure- 
dilatation 

Pvd'  =  - a2pVM t  +  a3pe,M(  (17) 

We  will  now  use  the  DNS  data  to  verify  the  functional  dependence  stated  in  Eq.  (17)  and 
also  to  calibrate  the  model  coefficients  at  and  <23. 

Because  the  production  V  =  0  in  decaying  isotropic  turbulence,  the  variation  of  the 
incompressible  pressure-dilatation  with  e,  can  be  verified  using  DNS  of  isotropic  turbulence. 
The  ratio  p/'d'/(pe,A/(2)  is  shown  as  a  function  of  non-dimensional  time  in  Fig.  7a.  This 
ratio  reaches  an  equilibrium  value  by  a  time  of  0.25,  substantiating  the  validity  of  the  second 
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term  in  Eq.  (17).  Based  on  the  DNS  value  of  the  equilibrium  ratio,  the  model  coefficient  q3 
in  Eq.  (17)  is  taken  to  be  0.2.  The  remaining  part  of  the  model  for  the  pressure-dilatation 
is  calibrated  against  simulations  of  homogeneous  shear  flow.  Fig.  7b  shows  that,  in  accord 
with  our  model,  the  rapid  part  of  the  pressure-dilatation  scales  as  ~pP\l} .  The  ratio  (p1' d1  — 
0.2 ( s.Mf)/ (pV M?)  reaches  an  approximate  equilibrium  value  of  —0.1,  suggesting  that  the 
model  coefficient  Oj  =  0.1. 

Finally,  after  using  the  arguments  of  section  2  to  neglect  pc ' d'  relative  to  p1' d' .  the  model 
lor  the  pressure-dilatation  becomes 

p'd'  —  —ttipVMf  +  OzptaMf  (IS) 

where  V  =  -  ff, ./(/,''(/"  is  the  production  of  kinetic  energy,  ,U<  =  \J'2l\ / ~;IU'  the  turbulent 
Mach  number.  es  the  solenoidal  dissipation,  and  the  model  coefficients  are  o2  =  0.1.  o3  =  0.2. 
We  are  now  in  t  he  process  of  applying  the  pressure-dilatat  ion  model  to  the  compressible  shear 
layer  and  the  flat  plate  boundary  layer. 

4  Conclusions 

We  have  obtained  a  model  for  the  pressure-dilatation  after  applying  scaling  arguments  ap¬ 
plied  to  a  formal  solution  for  homogeneous  turbulence.  This  model  has  been  validated 
by  direct  comparison  with  DNS  results  for  pressure-dilatation  in  homogeneous  shear  flow 
and  isotropic  turbulence.  Eq.  (18)  is  a  reasonable  approximation  for  inhomogeneous  flows 
without  walls.  However,  future  refinements  may  be  necessary  to  capture  different  physical 
processes  of  importance  in  other  flows.  For  example,  a  process  of  importance  in  shock- 
turbulence  interaction  and  engine  flows  is  a  compressive  mean  velocity  field.  In  such  flows, 
the  pressure-dilatation  model  Eq.  (18)  will  have  a  contribution  which  is  linear  in  the  mean 
compression.  Inspection  of  the  dilatation  equation  for  homogeneous  compression,  obtained 
by  taking  the  divergence  of  the  momentum  equation,  shows  that  there  should  be  such  a 
term  in  the  pressure-dilatation.  Of  course,  even  though  the  form  of  the  dependence  on  mean 
compression  is  already  present  in  our  model,  in  order  to  get  the  coefficient  of  the  dependence 
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right.  il  may  be  necessary  to  to  add  a  term  like  a  {J>K iim  mM*  to  the  model,  with  a  ,  perhaps 
being  calibrated  with  respect  to  DNS  of  homogeneous  compression.  In  wall-bounded  (lows, 
the  pressure-dilatation  is  probably  smaller  than  the  estimate  given  by  hq.  (IS)  because,  first, 
as  shown  by  Kim  and  Lee10  the  rms  rapid  pressure  (except  very  near  the  wall)  is  a  smaller 
fraction  of  /></*  in  wall  bounded  flows  compared  to  homogeneous  shear  <hnv.  and  second,  the 
velocity  component  normal  to  the  wall  is  preferentially  damped. 

Ihe  model  presented  Imre  lor  the  pressure-dilatation  is  currently  being  applied  in  the 
calculation  of  inhomogeneous  hows.  As  pointed  out  in  the  previous  paragraph,  some  exten¬ 
sions  may  be  necessary  in  ordei  to  incorporate  different  physical  processes  in  other  Hows. 
We  plan  to  further  develop  th  pressure-dilatation  model  for  more  general  (lows. 
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Table  i:  Parameters  for  the  DNS  cases  of  homogeneous  shear  how 


Case 

s 

V 

Mt.  o 

Ry  o 

(>' 

cl' 

SI 

10 

1/150 

0.3 

24 

0 

0 

S2 

15 

1/150 

0.3 

24 

0 

0 

S3 

15 

1/125 

0,1 

20 

0 

Table  2:  Parameters  for  the  DNS  cases  of  decaying  isotropic  turbulence 


( 'ase 

l > 

Mt.  o 

Ry  o 

h' 

d! 

1)1 

1  /2  00 

0.4 

ilh 

0 

0 

BB 

1/175 

0,5 

0 

6» 

1/175 

0.6 

27 

0 

0 

1 
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Fig.  2a.  Comparison  of  p  d!  and  p  d'  for  Case  SI. 


Fig.  2b.  Comparison  of  the  time  integrals  of  p~d'  and  pr  d'  for  Case  SI. 
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Fig.  3a.  Comparison  of  pa'  and  pr d!  for  Case  S2. 


/w  I* 

Fig.  3b.  Comparison  of  the  time  integrals  of  p  d!  and  p  d!  for  Case  S2. 
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Rg.  4a.  Comparison  of  p  d'  and  p  d!  for  Case  S3. 


Rg.  4b,  Comparison  of  the  time  integrals  of  p  d!  and  p  d!  for  Case  S3. 
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Fig.  7a.  Incompressible  pressure-dilatation  in  decaying  isotropic  turbulence. 
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Fig.  7b.  Incompressible  pressure— dilatation  in  homogeneous  shear  turbulence. 
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